Simultaneous Measurement of Ionization and Scintillation from Nuclear Recoils in 
Liquid Xenon as Target for a Dark Matter Experiment 
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We report the first measurements of the absolute ionization yield of nuclear recoils in liquid xenon, 
as a function of energy and electric- field. Independent experiments were carried out with two dual- 
phase time projection chamber prototypes, developed for the XENON Dark Matter project. We 
find that the charge yield increases with decreasing recoil energy, and exhibits only a weak field 
dependence. These results are a first demonstration of the capability of dual phase xenon detectors 
to discriminate between electron and nuclear recoils, a key requirement for a sensitive dark matter 
search at recoil energies down to 20 keV. 

PACS numbers: 95.30.Cq Elementary particle processes in Astrophysics , 95.35.+d Dark Matter , 25.40.Dn 
Elastic neutron scattering , 29.40.Mc Scintillator Detectors 



Introduction: Current evidence indicates that one 
quarter of the mass-energy density of the universe is 
composed of cold, non-baryonic dark matter, which has 
thus far been observed only through its gravitational 
interactions with normal matter. Its precise nature is 
undetermined, but weakly interacting massive particles 
(WIMPs) are an attractive candidate which may be de- 
tectable via rare elastic scattering interactions deposit- 
ing a few tens of keV in target nuclei. For a review 
of the motivation for WIMPs, and current experimental 
WIMP searches, see The best limits, from CDMS, are 
<0.06 events/kg/day (in Ge) 2]. Improving the search 
sensitivity will require both larger detectors and lower ra- 
dioactive backgrounds. An attractive method to achieve 
this uses liquid noble gas-based detectors, which promise 
to be readily scalable to the multi-ton scale. For the 
XENON experiment |3( , we are developing a large volume 
liquid xenon (LXe) dual-phase time-projection-chamber 
(TPC) that features simultaneous measurement of recoil 
ionization and scintillation to determine the energy and 
3-D localization on an event-by-event basis. Nuclear re- 
coils from WIMPs (and neutrons) have denser tracks, 
and thus have been assumed to have greater electron-ion 
recombination than electron recoils, providing a basis for 
discrimination against radioactive background gammas 
and betas. However, ionization from these nuclear re- 
coils has not been measured until now. In this report, we 
describe the first measurements of the ionization yield of 
low energy nuclear recoils in LXe and confirm the base- 
line discrimination capability of this technique. 

Independent measurements were carried out with two 
separate dual-phase detectors at Case and Columbia. 
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The active LXe is defined by two electrodes, a cathode 
and a gate grid, at a distance of 1.9 cm (Columbia) and 
1.0 cm (Case). An event in the active LXe produces 
both ionization electrons and scintillation photons. If 
trapping by impurities is negligible, the electrons which 
escape initial recombination drift freely under an applied 
electric field of a few kV/cm up to the liquid-gas phase 
interface where they are extracted and accelerated by 
the high (~10 kV/cm) field applied across the gas gap. 
These accelerated electrons produce proportional scintil- 
lation photons before being collected on the anode. The 
fields are created by setting potentials on four electrodes: 
cathode, gate grid, anode, and top grid. The electrodes 
are stretched wire grids, except for a solid stainless steel 
cathode in the Case detector and a mesh top electrode 
in the Columbia detector. The wire grids are made with 
120 fim diameter BeCu wires on a 2 mm pitch, except for 
40 fim gate grid in the Case detector, and are soldered 
to SS (Columbia) and Cirlex (Case) frames. 

Photomultiplier tubes (PMTs) detect both the direct 
scintillation light (SI) produced at the event site in the 
liquid and the proportional scintillation light (S2) pro- 
duced in the gas. The time difference between the two 
signals, determines the depth of the event. The 5 cm 
diameter, 4 cm long metal channel PMTs (Hamamatsu 
Model R9288) were developed to be directly coupled with 
LXe, and have a quantum efficiency of ~15% at the 178 
nm wavelength of the Xe scintillation 4]. PTFE is used 
as a VUV reflector enclosing the active LXe volume in 
both detectors. The Columbia detector uses two PMTs, 
one in gas and one in liquid. Because of total internal 
reflection at the liquid surface (index of refraction 1.65 
at 187 K 5]), the light collection is significantly better 
for the PMT in the liquid. The total SI light collection 
efficiency (fraction of photons incident on the photocath- 
ode of a PMT) is estimated from Monte Carlo to be more 
than 50%, uniform across the volume to within 5%. An 
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internal blue light LED is used to calibrate and monitor 
the gain of the two PMTs. The Case detector uses a 
single PMT in the gas, with an SI light collection effi- 
ciency of ~16%. The PMT gain is monitored using the 
SI light from 5.3 MeV alphas from a source described 
below. The operation of a dual phase LXe detector is 
described in more detail in |6|. 

Both detectors use a liquid nitrogen-cooled cold finger 
cryostat with temperature control to maintain a stable 
liquid temperature. The detectors were operated with 
vapor pressures between 2.0 and 2.8 atm, with better 
than 1% stability. The Case detector also uses a triple- 
parallel plate capacitor system to align the liquid surface 
with the grids and monitor the level and stability of the 
liquid surface, which waspj stable to within 20 fim over 
two months. The LXe in both detectors was purified us- 
ing a high temperature getter. The Columbia detector 
used the purification system with continuous gas circula- 
tion through the getter, developed for the first XENON 
prototype |3, Q • The Case detector used a similar recir- 
culating system only at the start of its two-month run. 

External gamma ray sources, including 57 Co, 133 Ba, 
and 137 Cs, were used for calibration of both detectors. 
The Case detector also had an internal 210 Po source de- 
posited on the center of the cathode plate, providing 
5.3 MeV alpha particles and 100 keV 206 Pb nuclear re- 
coils. Events <0.5 mm above the cathode were tagged 
as possible alphas and 206 Pb recoils. For neutron data, a 
5 Ci AmBe source (Columbia) and 25 /iCi 252 Cf source 
(Case) were used, with lead shielding to attenuate gam- 
mas from the sources. 

The PMT signals were digitized with multiple ADCs 
sampling at 5 MHz-1 GHz. The effective trigger thresh- 
old for the Case detector was 4.5 electrons for the S2 sig- 
nal, and 50% single photoelectron (spe) acceptance for 
SI. For the Columbia detector the trigger was either the 
coincidence of the SI signals from the two PMTs (at a 
few spe level) or the SI signal from the PMT in the liquid 
(at a ~6 spe level). 

Calibration and Raw Data: Both detectors were 
calibrated primarily with 57 Co 122 keV gammas, for both 
the SI and S2 signals. The SI signal serves as an energy 
axis for nuclear recoils, given by E r = E e /C e ff • S e /S r 
where E r is the nuclear recoil energy, E e is a linear elec- 
tron recoil energy scale based on the 57 Co SI peak, C e ff 
is the effective Lindhard factor relating the scintillation 
yields of nuclear and electron recoils at zero field, and S x 
is the loss of scintillation light due to the recombination 
suppression by the electric field for species x, equal to 
S(E)/So in Fig. 0] Recent measurements at Columbia 
give S r for 56.5 keV recoils and C e jfdX e nerg ies up to 
56.5 keV 10]. Other experiments [ill E El El have 
measured C e ff at higher energies. The Columbia and 
other data, except 0, were parameterized by C e ff = 
0.0984£ r - 169 . For simplicity, we assume S r has no en- 
ergy dependence since the Columbia measured value is 
close to unity, as seen in Fig. 21 Also shown in Fig. 21 
are measurements of S e as a function of drift field. Here- 



after, nuclear recoil energies calculated in this way will 
be denoted 'keVr', and electron recoil energies based on 
the linear SI scale will be denoted 'keVee'. In the Case 
and Columbia detectors the SI calibration gives, respec- 
tively, 1.5 photoelectrons (pe)/keVee and 5 pe/keVee, at 
zero electric field, which correspond to 0.28 pe/keVr and 
1 pe/keVr for 55 keVr nuclear recoils. 

The S2 signal measures the ionization from each event. 
There is no published data of charge yield vs field for low 
energy gammas in LXe, so this measurement was per- 
formed with 57 Co 122 keV gammas, in both the Case 
[TEj and Columbia detectors operated in single (liq- 
uid) phase. The results are shown in Fig. 0] The S2 
calibration is a function of gas pressure and electric field 
across the gap [l?|, with the added complication of pos- 
sible electron multiplication near the anode wires. Typ- 
ical values were 19 pe/e - with 4.6 kV/cm (Case) and 
8.4 pe/e" with 2 kV/cm (Columbia). 

Figures Hand [21 shows the detectors' responses to neu- 
tron and low energy Compton scattering events. The 
logarithm of the ratio S2/S1 is plotted as a function of nu- 
clear recoil energy (keVr). In both detectors, the elastic 
nuclear recoil band is clearly separated from the electron 
recoil events. The events identified around 40 keVee are 
from gamma rays produced by neutron inelastic scatter- 
ing on 129 Xe. Additional gamma ray events are emitted 
following inelastic scattering of neutrons on 131 Xe and 
19 F contained in PTFE. The detector responses to 137 Cs 
(Columbia) and 133 Ba (Case) show only the gamma ray 
band. 

Nuclear Recoil Ionization Yield: The ionization 
yield is defined as the number of observed electrons per 
unit recoil energy (e~/keVr). The nuclear recoil ioniza- 
tion yield as a function of energy is shown in Fig. 
for several drift fields. The uncertainty on the yield is 
dominated by the systematic error from the 57 Co S2 cal- 
ibration. 

Fig. Q] summarizes the relative light and charge yields 
as a function of drift field for different particles in LXe: 
122 keV gamma rays from [HJ 56 keV Xe nuclear 
recoils from wl and this paper, and 5.5 & 5.3 MeV al- 
phas from [3 an d the Case detector. The relative light 
yield S(E)/So is the light yield relative to that at zero 
field So. The relative charge yield Q(E)/Qq is the charge 
collected relative to that at infinite field (i.e., with no re- 
combination) Qq. For nuclear recoils, Qo = E r • C/W e , 
where E r is the recoil energy, C is the supression pre- 
dicted by Lindhard theory |l<|, and W e = 15.6 eV Hj 
is the average energy required to produce an electron-ion 
pair in LXe. For gammas and alpha particles, we divide 
the known energy by W e to obtain Qo- 

Discussion: The important characteristics of the 
measured ionization yield of nuclear recoils in LXe are 
its value relative to that of other particles (Fig. 0}, its 
energy dependence (Fig. |3J), and its field dependence 
(Fig. 0J. The yield is a function of the Lindhard fac- 
tor and the amount of recombination in the track. The 
Lindhard factor has no field dependence, and does not 
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FIG. 1: Columbia detector response to AmBe neutron (top) 
and 137 Cs gamma sources (bottom), at 2 kV/cm drift field. 
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FIG. 3: Energy dependence of nuclear recoil ionization yield 
at different drift fields. 
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FIG. 2: Case detector response to 252 Cf neutron (top) and 
133 Ba gamma sources (bottom) at 1.0 kV/cm drift field. 



FIG. 4: Top: Field dependence of scintillation and ionization 
yield in liquid xenon for 122 keV electron recoils (ER), 56 
keVr nuclear recoils (NR) and alphas. Bottom: Ionization 
yields scaled by their 4.5 kV/cm values. 

affect the relative charge yield (Fig. Q}. Lindhard does 
predict a slight decrease in charge yield with decreasing 
energy, the opposite of what we see. Therefore we expect 
recombination and its dependence on track density and 
shape to explain the above phenomena. 

Recombination is primarily a function of electric field 
and ionization density, with stronger recombination at 
low fields and in denser tracks. Ionization density along 
a track corresponds roughly to electronic stopping power, 
plotted in Fig. [51 for alphas, electrons, and Xe nuclei in 
LXe, as given by ASTAR, ESTAR, and SRIM [111, re- 
spectively . Also shown is the total energy lost to elec- 
tronic excitation per path-length for Xe nuclei, which 
differs from the electronic stopping power in that it in- 
cludes energy lost via electronic stopping of secondary 
recoils |22j. 

At lower energies, the stopping power for Xe nuclei de- 
creases, indicating lower recombination and higher ion- 
ization yield, as observed. The stopping power, including 
daughter recoils for Xe nuclei at 56 keV is higher than 
that of alphas at 5.5 MeV. This is in conflict with what 
we observe in the relative ionization yields, indicating 
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FIG. 5: Predicted electronic stopping power, dE/dx, for dif- 
ferent particles in LXe. The circles indicate the particle en- 
ergies discussed in the text. 



that the higher ionization yield of a low energy Xe re- 
coil compared to an alpha is due to the different track 
geometries for these two particles. Recombination for al- 
phas is enhanced because of their cylindrical track shape 
with a very dense core and a "penumbra" of delta rays. 
SRIM simulations show that nuclear recoil tracks have 
many branches, with most energy lost by secondary re- 
coils. Each of these secondary branches presumably end 
in a very sparse track since the stopping power falls with 
energy. This may result in a low-density halo of ioniza- 
tion around the track, which does not recombine even 
at very low fields. If this lack of recombination extends 
to zero field, it would explain the difference between the 
measured scintillation yield for nuclear recoils and the 
Lindhard prediction. The importance of the track ge- 
ometries is also seen in the striking contrast between the 
field dependence of nuclear recoils and those of electrons 
and alphas (bottom graph of Fig. QJ, though electrons 
have a much lower, and alphas a comparable stopping 
power to nuclear recoils. Our new finding of the weak 
field dependence for nuclear recoils in LXe has the im- 
portant practical consequence for dark matter detection 



that only modest fields are needed for the background- 
discriminating measurement of nuclear recoils. 

Background Discrimination: The different values 
of S2/S1 between nuclear recoils and electron recoils 
make liquid xenon an excellent target to discriminate be- 
tween these two types of interaction. Acceptance win- 
dows for nuclear recoils were defined from the neutron 
elastic recoil band as shown in Fig. ^ an d [21 The re- 
jection of background electron recoils was measured by 
irradiating both detectors with pure gamma ray sources 
( 133 Ba or 137 Cs), which produce low energy electron re- 
coils from Compton scatters, and recording the number 
of electron recoils outside of 50% acceptance nuclear re- 
coil windows. The Columbia detector, with higher SI 
light collection, achieved an electron recoil rejection effi- 
ciency of 98.5% independent of energy down to 20 keVr 
at 2 kV/cm drift field. This efficiency appears to be lim- 
ited by charge loss at the edge of the detector where the 
field was non-uniform (no field shaping rings were used 
in this prototype). The rejection efficiency in the Case 
detector improved with increasing energy (because of the 
non-optimal light collection) and was better than 99.5% 
above 80keVr, presumably due to a better field geometry 
with closer electrodes and better S2 resolution. The re- 
jection efficiency improves with higher drift field (in the 
Case detector, e.g., from 91% at 0.2 kV/cm to 97% at 4.6 
kV/cm for 35 keVr). With an optimized field configura- 
tion, and with the XY position sensitivity of a LXeTPC 
as proposed for the XENON experiment 3], the prob- 
lem of edge events should be negligible. The dual phase 
xenon technique thus appears very promising for a large- 
scale detector with powerful electron recoil background 
discrimination for a sensitive WIMP dark matter search. 
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